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Strongly-coupled gauge theories are an important ingredient in the construction of many exten- 
sions of the standard model, particularly for models of electroweak symmetry breaking in which 
the Higgs boson is a composite object. There is a large parameter space of such gauge theories 
with a rich phase structure, particularly in the presence of large numbers of fermionic degrees of 
freedom. Lattice simulation provides a non-perturbative way to explore this space of strongly- 
coupled theories, and to search for interesting dynamical features. Here I review recent progress 
in the simulation of such theories, with an eye towards applications to dynamical electroweak 
symmetry breaking. 
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1. Introduction 

Lattice gauge theory has enjoyed considerable success as a non-perturbative method to study 
and understand quantum chromodynamics (QCD), the only example of a strongly coupled gauge 
theory observed in nature thus far. With modern computers and computational techniques, lattice 
QCD is attaining high levels of precision, with simulation parameters approaching the physical 
point. However, the success enjoyed by lattice QCD is no reason for us to restrict the attention of 
our lattice studies to that theory alone. QCD is far from unique as an example of an asymptotically 
free Yang-Mills gauge theory, inhabiting a large space of models, all of which are well-suited to 
numerical study on the lattice. Some of these models have properties very different from QCD, 
lacking phenomena such as confinement and spontaneous breaking of chiral symmetry, and instead 
showing the restoration of conformal symmetry at low energies. 

These more general Yang-Mills theories are interesting in their own right, particularly in the 
context of understanding the implications of conformal symmetry in four-dimensional quantum 
field theory and the AdS/CFT correspondence. But in the era of the Large Hadron Collider, our 
primary motivation for the study of theories other than QCD is provided by the search for physics 
beyond the standard model (BSM). Strongly-coupled field theories are an essential ingredient for 
models of dynamical electroweak symmetry breaking, triggered by the spontaneous breaking of 
chiral symmetry by the new gauge sector. This can be accomplished either directly with no scalar 
Higgs field required (as in technicolor [1, 2]), or by generating a composite Higgs as a pseudo- 
Goldstone boson (as in Little Higgs theories [3, 4]). I will focus on technicolor models as a moti- 
vation in this review, since the presence of strongly-coupled physics close to the electroweak scale 
makes non-perturbative input for such models more urgent. 

Such strongly-coupled models can be difficult to work with, as many quantities cannot be 
calculated using perturbation theory. Early work on technicolor relied heavily on QCD experiment 
for non-perturbative input, leading to strong tensions with experiment through the S'-parameter 
[5] and flavor-changing neutral cun^ents [6, 7]. But there is no reason to expect more general 
strongly-coupled theories to lead to the same problems. Indeed, for many-fermion gauge theories, 
it has been hypothesized that a theory with approximate scale invariance (colloquially known as a 
"walking" theory, since the running gauge coupling evolves relatively slowly) should have the right 
properties to cure both the flavor-changing neutral current problem [8, 9, 10] and the S'-parameter 
[11]. Lattice simulation provides us with an ideal method to confirm or refute these conjectures 
non-perturbatively, and to search for other interesting dynamical changes as we move away from 
QCD and towards the unknown. (For a review of lattice BSM work with greater focus on these 
phenomenological issues, see [12].) 

The primary focus of this review, and the subject of a recent surge in interest from a number 
of lattice groups, is on the phase structure of Yang-Mills gauge theories with large numbers of 
massless fermions. In section 2, I introduce the concept of the "conformal window," a range of 
theories which approach an infrared fixed point at low scales. Some limited insights from perturba- 
tion theory and other analytic arguments are presented. Section 3 reviews various lattice results for 
"QCD-like" theories that are known to lie below the conformal transition. In section 4, techniques 
and simulation results for theories which may be inside the conformal window are discussed. Fi- 
nally, in section 5 I provide some concluding remarks. 
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2. The conformal window 

QCD is our most familiar example of a strongly-coupled Yang-Mills gauge theory, but it is 
far from unique. Even restricting our attention to SU(A'^c) gauge symmetry groups, there is a vast 
parameter space available by varying the number of colors Nc and the fermion content of the theory. 
Our starting point is the Yang -Mills Lagrangian, which for a theory with Nf massless Dirac fermion 
species in representation R of the gauge group reads 

Nf 

= -TiV^" +"Lwmwi, (2.1) 

where D is the usual gauge-covariant derivative. More generally one may include fermions in 
multiple representations of the gauge group [13], but I will consider only the single-representation 
case here. 

At the classical level this Lagrangian is scale invariant, but quantum effects break the confor- 
mal symmetry, with the resulting dependence on renormalization energy scale pt appearing in the 
renormalized coupling constant g{pL). The scale dependence of g{pL) can be computed order-by- 
order in perturbation theory, and is usually expressed in terms of the jS -function: 

= 37r^ = -i8o«' - - i82«' - ... (2.2) 
with a{pL) = g{^Y /An. The universal values for the first two coefficients are 



34 / 20 \ 

-N^, - [aC2{R) + yA^.J T{R)Nf 



(2.4) 



with T{R) and C2{R) the trace normalization (first Casimir invariant) and quadratic Casimir invari- 
ant of representation R, respectively. A list of invariants for some representations of SU(A'^) are 
shown in table 1; a more exhaustive tabulation of group-theory factors can be found in [14]. This 
expansion can of course be continued to higher order, at the price of specifying a renormalization 
scheme, such as the MS scheme in which the first four coefficients are known [15]. 

We can observe several features in the Nc-Nf plane based on the two-loop j8 -function given 
above, for some fixed representation R. For any given value of A'^ , there is a critical number of 
flavors 

Nf = (2.5) 

at which the first coefficient j3o changes sign. Theories with Nf > Nf^ are no longer asymptotically 
free (the trivial solution j3 (0) =0 is no longer ultraviolet-attractive), and will not be considered 
further here. If instead Nf is taken to be lai^ge but slightly below Nf'^, then for any representation 
there is a range of values satisfying j8i < but jSq > 0. The two-loop jS -function then has a second 
fixed-point solution [16, 17] at non-zero coupling strength: 

af'^ = -§'. (2.6) 
Pi 
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Table 1: Casimir invariants and dimensions of some common representations of SU{N): fundamental (F), 
two-index symmetric (52), two-index antisymmetric (A2), and adjoint (G). 



So long as the coupling a; is sufficiently weak, the presence of this infrared-stable fixed point 
ensures that the theory will be perturbative at all scales, so that the use of the two-loop j3 -function 
is self-consistent. 

As Nf is decreased from Nj^, the perturbative fixed-point coupling increases, and the use 
of perturbation theory ceases to be valid. Still, we know that the non-perturbative properties of 
confinement and spontaneous chiral symmetry breaking should be recovered at small enough Nf, 
in the regime of QCD {Nf = 2) and the quenched {Nf = 0) limit. We define a critical value Nj- as 
the point where the infrared behavior of the theory changes from conformal to confining, so that 
the range of theories 

N'f < Nf < — — ^ (2.7) 

all possess infrared fixed points. This range is known as the "conformal window". Perturbation 
theory is almost certainly unreliable to describe physics at the infrared fixed point once Nf is close 
to Nj, so that study of this transition requires a non-perturbative approach. A sketch of the behavior 
outlined here in the Nc-Nf plane is shown for the fundamental representation in fig. 1. 

Due to the presence of an infrared fixed point, theories inside the conformal window are not 
expected to exhibit confinement or spontaneous chiral symmetry breaking. Conversely, no theory 
with spontaneously broken chiral symmetry can be in the conformal window, since the infrared 
fixed point is unstable with respect to the fermion mass m, a relevant perturbation. Another way to 
state this point is that if the fermions develop a mass m 7^ 0, then in the infrared limit they will be 
screened out of the theory [18], and the coupling will run as in a pure-gauge theory at low scales. 

As I have stressed, although perturbation theory is useful for introducing the basic concepts of 
the conformal window, it ceases to be a useful description near the lower edge. Lattice simulation 
provides a fully non-perturbative, albeit numerical, approach to the study of these near-critical the- 
ories, and can allow us to constrain Nf. Early pioneering studies [19, 20, 21] were entangled with 
lattice-artifact phase transitions and unable to draw clear conclusions, but a growing number of lat- 
tice groups have undertaken many-fermion simulations over the past few years, and our knowledge 
of the conformal window in certain classes of theories is rapidly improving. 

In fig. 2, 1 summarize the efforts of lattice groups thus far in attempting to study many-fermion 
theories. Each box on the plot represents a particular theory, i.e. a particular choice of the pa- 
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Figure 1: Cartoon showing the phase structure in the Nc~Nf plane for SU(A^c) Yang-Mills with fundamental 
fermions. The region labelled "AF lost" contains theories with /3o < 0, which are not asymptotically free. 
"CBZ" denotes Caswell-Banks-Zaks theories [16, 17], which flow to very weak infrared fixed points. The 
edge of the shaded "conformal" region, set by the condition al < 1, is intentionally faded to indicate that 
perturbation theory is becoming unreliable. Non-perturbative methods must be used to find the true edge of 
the conformal window. 



rameters {Nf,Nc,R), and each point represents the combined efforts of one lattice group for that 
theory. Simulation results are color-coded based on whether their results indicate confinement 
(red), conformality (blue), or are inconclusive (puiple), based on my reading of the authors' own 
claims. Theories themselves are colored similarly, based on my own (subjective) determination 
of the overall consensus in the field. Most results thus far are focused on the Nc = 3,R = F case, 
which contains QCD; these theories allow for the reuse of codes written for QCD with minimal 
modification. 

I have included two analytic approaches to determining Nj on this plot. A direct estimate, 
denoted by a dashed green line, is obtained by studying the Schwinger-Dyson equation for the 
fermion propagator [22]. To make progress in this approach, several assumptions must be made, 
including a truncation of the set of diagrams contributing to the fermion self-energy known as 
the "ladder approximation". Although some qualitative understanding of the physics behind the 
conformal transition may be gained with this approach, the error introduced by this truncation is 
not well controlled at strong coupling, and the perturbative estimates of Nj require the further use 
of the two-loop j3 -function. They are included only for reference. 

The other analytic result indicated in this summary plot is a conjectured bound on Nj known 
as the "thermal inequality" [23]. The underlying assumption of this inequality is that the number 
of degrees of freedom should not increase as a theory flows from ultraviolet to infrared; this is an 
intuitive statement if one considers the RG flow in terms of a "blocking" operation which tends to 
reduce the complexity of the system. The thermal inequality is based on the free energy, which in 
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the free-field limit of a theory with scalars, A^v vectors and Nf Dirac fermions is equal to 



27-4 



FiT) 



"90" 



A^, + 2A^, + -(4A^j) 



It A 



~90 



■/(r). (2.8) 



The function f{T) gives a count of the number of light degrees of freedom, so that the conjecture 
reads 

lim/(r) = //«</c/i/= lim/(r). (2.9) 



This inequality can be used to bound N'j by applying to a theory in the chirally broken phase. In 
the infrared limit, the relevant degrees of freedom for such a theory are the Nj — l Goldstone bosons 
associated with the breaking pattern S\J{Nf) x S\J{Nf) — > S\J{Nf). Comparing to the ultraviolet 
count of free fermions and gauge bosons, we arrive at the result 



/ 16 

implying that A'^ < 4Nc. It should be noted that this relation applies only for Nc > 3; Nc = 2 is a. 
special case. Because the representations of SU(2) are real or pseudo-real, the bilinear (yT^W) 
is invariant under an enhanced chiral symmetry group S\J{2Nf), rather than the usual SU(A^/)x 
SU(A'y ). The breaking pattern preserving the most symmetry is then S\J{2Nf) — )• Sp(2A'y ), chang- 
ing the counting for the IR degrees of freedom. The resulting bound for Nc = 2is then N'^ < 4.7. 

In two dimensions, the conjecture eq. (2.9) is proven for an RG flow connecting two conformal 
field theories, using Zamolodchikov's c-theorem [24] and the direct relation of central charge c to 
the free energy. A proof for the generalization of the c-theorem to d = 4 known as the a-theorem 
has recently been presented [25], but its connection to the free energy is less straightforward, so 
an a-theorem proof does not imply the thermal inequality. A distinct constraint based on the trace 
of the energy-momentum tensor in the UV and IR limits does follow from the a-theorem, but the 
resulting bound on Nj^ is very weak, lying above A^^^ [26]. 



3. QCD-like theories 

In using the term "QCD-like", I am referring specifically to theories with far-infrared proper- 
ties resembling QCD, namely confinement and spontaneous breaking of chiral symmetry. These 
are desirable properties in the context of technicolor or composite Higgs theories, as the broken 
chiral symmetry is generally used to trigger electroweak symmetry breaking. (It is also possible to 
start with a theory inside the conformal window and give masses to some of the fermions, triggering 
spontaneous symmetry breaking for the massless species, as in e.g. [27, 28].) 

Because their infrared properties are similar, lattice simulations of QCD-like theories can rely 
on mature tools and techniques developed for use in QCD itself. In particular, the presence of chiral 
symmetry breaking means that the chiral Lagrangian provides a good effective low-energy descrip- 
tion of any QCD-like theory. Chiral perturbation theory is an essential framework for taking lattice 
results at the relatively heavy quark/fermion masses which are often required for good control of 
finite-volume effects, and extrapolating to reach the physical point (for QCD) or the massless limit 
(in our current context). 
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Figure 2: Summary plot of recent lattice simulations searching for the conformal window and its properties 
for various choices of A^^. and fermion representation. Each box represents a particular theory {Nc,R,Nf), 
and each dot corresponds to the combined results of one lattice group. 
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SU(2) and SU(3) chiral perturbation theory are familiar to most lattice theorists, and the gen- 
eralization to arbitrary SU(A'y ) is straightforward [29], and has been carried out to next-to-next-to- 
leading order for several observables [30, 31, 32]. Because the higher-order terms contain explicit 
loop contributions linear in Nf, the convergence of chiral perturbation theory is worse with more 
fermions; the expansion pai^ameter of the theory scales as NfMp/Fp (up to numerical factors.) 
For a fixed lattice spacing a, then, we expect that the bare fermion mass should be decreased as 
m oc Y^A^ in order to maintain the stability of the chiral expansion. 

The set of theories with Nc = 3 and fundamental-representation fermions is a natural choice 
for the exploration of QCD-like theories, because it contains QCD itself as a reference point. Fur- 
thermore, this choice allows the re-use of existing lattice codes which are heavily optimized for 
manipulation of SU(3) matrices. I will focus in this section on theories from Nf = 2 to Nf = 8, 
which by consensus of existing studies are outside the conformal window (as shown in fig. 2). 

A relatively conservative approach has been taken by the LSD collaboration, who have gen- 
erated a set of Nf = 6 and Nf = 2 gauge ensembles, tuned so that confinement-scale observables 
(ro, nip, etc.) are matched in lattice units between the two sets. They have computed a number 
of different observables which are particularly interesting in the context of technicolor theories, 
including the chiral condensate [33] and 5-parameter [34]. 

Some additional details of the latter calculation, which saw a reduction in S for the Nf = 6 
theory compared to naive scaling predictions, have been released in these proceedings [35]. The 
issue of finite- volume corrections, a possibly significant source of systematic error in the previous 
results, is investigated in detail. By studying the same measurement on smaller volumes (L/a = 
16), an artificial reduction in S due to finite-volume effects is accompanied with a "freezing" of 
the pseudoscalar mass which is not observed on the larger {L/a = 32) lattices. Furthermore, a 
"pole-dominance" estimate of 5 which uses sum rules and the properties of the lowest-lying vector 
and axial- vector states shows a spurious reduction on the larger L/a = 32 volume for the Nf = 2 
estimate, whereas no reduction is seen in the direct computation of S. These results imply that S 
may be less sensitive to finite-volume effects than other quantities, and that the reduction seen at 
Nf = 6 is physical. 

Another quantity recently computed on the Nf = 2 and Nf = 6 lattices is the (A/ = 2) n-n 
scattering length [36]. In the context of a technicolor theory, where the longitudinal degrees of the 
W and Z bosons are technipions, this calculation is related to WW scattering, a process which is 
experimentally accessible at the LHC. Results for the low-energy constants a4 and as in the Nf = 2 
theory gives values which are likely too small to be probed even with 100 fb^^ of LHC data [37]. 
Extending the analysis to Nf = 6 faces further complications, as the appearance of additional low- 
energy constants in the scattering length makes it impossible to isolate the parameters of interest. 
However, an increase in the overall next-to-leading order contribution to the scattering length is 
observed from Nf = 2 to Nf = 6. If this qualitative trend continues as Nf increases, then the effect 
may become observable in the full LHC data set. Further study will be needed, not only at different 
Nf but also in different tt-tt scattering channels. 

Moving towards the conformal transition, we turn now to the Nf = S theory. Several inde- 
pendent groups using different methods have shown that Nf = 8 is outside the window [38, 39, 
40, 41, 42, 43, 44, 45, 46, 47, 48], so that techniques appropriate for QCD-like theories should 
still be applicable. Furthermore, Nf = 8 has a direct counterpart in the technicolor literature, the 
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"one-family" model [1] (although because the technifermions in this model cany ordinary color, 
it is severely constrained by existing LHC searches [49]). It will be quite interesting to see if the 
trends indicated atNf = 6 continue as the transition is approached. 

However, few simulations have been carried out thus far with a focus on phenomenologically 
relevant observables in this theory. The authors of [45] have carried out simulations of the Gold- 
stone sector at Nf = S, but are not able to obtain consistent fits to chiral perturbation theory with 
the available range of masses and volumes. The LatKMI collaboration has recently begun simula- 
tions with the highly improved staggered quark (HISQ) action at every Nf multiple of 4, including 
Nf = 8 [50]. Their results, again focusing on the Goldstone sector, are too preliminary to draw any 
strong conclusions, but the authors report a clear qualitative difference between Nf = 4 and Nf = 8. 

Another group has attempted to extract the deconfinement temperature Tc [47, 48]. They ob- 
serve a finite-temperature transition at Nf = 0, 4, 6, and 8, based on results for the chiral condensate 
and Polyakov loop. For the Nf = 6,8 theories, simulations at different A'; are carried out, in order 
to show the expected scaling behavior with N for a real deconfinement transition, T^. ~ \/a{^c)Nt- 
(Theories within the conformal window also show deconfinement transitions on the lattice, but they 
are lattice artifacts and do not persist in the continuum limit; see e.g. [51].) The authors go on to 
compute Tc/Kref, where A^/ is a "UV reference scale" defined through the two-loop j8 -function. 
The scale A^/ is dependent on the choice of a reference coupling strength, but for various such 
choices, the authors observe a decreasing ratio Td K-ef as Nf increases. This result is broadly con- 
sistent with analytic estimates of the scaling of Tc [52]. The question of the relationship of Tc to 
other more physical observables, e.g. bound-state masses, remains open. 

Finally, I turn to the case of A'^^ = 2 gauge theories (again with fundamental fermions), which 
have begun to attract more interest recently. Because the representations of SU(2) are real or 
pseudo-real, the bilinear (yj^y) invariant under an enhanced chiral symmetry group SU(2A^y), 
rather than the usual S\J{Nf)x S\J{Nf). A first study of the spectrum of Goldstone bosons in 
the Nc = 2, Nf = 2 theory has been carried out by Lewis, Pica and Sannino [53]. They show the 
expected number of five Goldstone modes, three pseudoscalar and two scalar, for the symmetry 
breakdown SU(4) — >-Sp(4), and additional bound states which remain heavy in the chiral limit. 
SU(2) theories hold paiticulai^ promise for the application of lattice methods, since the addition 
of four-fermion interactions and non-zero chemical potential is possible without giving imaginary 
contributions to the action. 

4. Infrared-conformal theories 

Although substantial progress can be made by studying the properties of QCD-like theories 
directly, eventually we must turn to the fundamental questions about the conformal window and 
the nature of the transition at its edge. Working with theories which are possibly governed by an 
infrared fixed point introduces many new challenges. Lattice simulations require the introduction 
of ultraviolet (lattice spacing) and infrai^ed (box size) cutoff scales, and typically include mass 
terms for the fermions as well. Great care must be taken to remove these introduced scales from 
the problem through continuum, infinite-volume and chiral extrapolations, in order to recover the 
approximate scale-invariance of the underlying theory. In some cases, controlled variation of these 
scales can be used to uncover certain properties of the continuum theory; we will discuss one such 



9 



Exploring Models for New Physics on the Lattice 



Ethan T. Neil 



example below, using bound-state masses in order to determine the mass anomalous dimension at 
an infrared fixed point. 

As I have done here, the concept of the conformal window is generally first introduced using 
the perturbative running coupling as a guide. A natural approach to studying such theories on the 
lattice is thus through non-perturbative definitions of a running coupling constant, or more broadly, 
studying the renormalization-group (RG) flow as a function of scale. This can be a daunting task in 
theories with many fermions, as the jS -function becomes numerically small, sometimes requiring 
evolution of the scale over many orders of magnitude in order to see significant RG flow. All 
such studies rely on the use of matching techniques to stitch together results obtained on lattices of 
different physical sizes (it is certainly impossible to change the scale by orders of magnitude while 
simulating at a fixed UV cutoff.) 

The Nc = 3, Nf = 16 theory provides a good starting point as a theory which is definitely within 
the conformal window. One of the earliest lattice calculations directed towards large-A^/ theories is 
[5 1], which studied the dependence of this theory on the bare gauge coupling j3 = 6/gQ. Despite the 
fact that A'/ = 16 is known to have a perturbative infrared fixed point, the authors see clear evidence 
for a phase transition as j8 is varied, with confinement and chiral symmetry breaking on the strong- 
coupling side, as indicated by lai^ge jumps in the string tension and chiral condensate, respectively. 
However, this transition is observed to occur at a fixed jSc independent of the simulation volume, 
indicating that it is associated with a lattice artifact and will not persist in the continuum limit. 

In a complementary study [54], Heller measured the Schrodinger Functional (SF) running 
coupling in the same theory at bare couplings somewhat below the transition. The coupling (L) 
was observed to decrease with increasing box size L, indicating that the j3 -function has the opposite 
sign relative to QCD. The simulations at these bare couplings were therefore on the strong-coupling 
side of the expected IR fixed-point, and therefore disconnected from the continuum limit at zero 
coupling. This observation of backwards RG flow is crucial in establishing the existence of a zero 
of the j3 -function using numerical methods with limited precision. 

Several different methods for determining the RG flow on the lattice, most of which focus on 
the extraction of a running coupling constant in a particular scheme, are now in use for a wide range 
of theories. 1 will not go into the details of these various methods here, but refer the interested reader 
to [55] for a thorough review. All of these techniques seek to establish the presence or absence of a 
zero in some analogue of the j3-function, with the observation of a sign change (i.e. reversal of RG 
flow) providing the clearest such signal. The approach of searching for a physical deconfinement 
temperature also falls into the class of RG-flow methods; the tell-tale signal of a theory outside 
the conformal window is given by the observation that Tc scales towards zero bare coupling as 
the temporal extent Nt is increased. As noted, an example of this method has been applied to the 
Nc = 3, Nf = S theory, indicating that it lies outside the conformal window [41, 47]. A list of 
references for all of the various RG-flow simulation results is given in table 2, and the resulting 
constraints on Nj are summarized in fig. 2. 

Another theory which has attracted considerable interest from lattice groups, and which is 
now generally agreed to be infrared conformal, is the Nc = 2 theory with Nf = 2 fermions in 
the adjoint representation. Simulations of this theory were first carried out in [56] with Wilson 
fermions, studying the pion decay constant and bound-state masses, but interpretation of the results 
were complicated by the presence of first-order bulk phase transitions, studied further in [57, 58]. 
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However, the authors did see some early hints that these scales were vanishing in the massless limit, 
interpreted as evidence for an infrared-conformal theory. 

In the context of this early work, it was not entirely clear how to extrapolate to the massless 
hmit. For a theory which is inside the conformal window, there is no spontaneous breaking of chiral 
symmetry, only explicit breaking introduced by the presence of a non-zero fermion mass m. At 
fixed m this breaking, which also leads to confining behavior in the infrared as the fermions screen 
out of the theory, leads to a spectrum which resembles a QCD-like theory, with the usual bound 
states. However, chiral perturbation theory will not apply in the description of the m-dependence 
of observables, due to the lack of spontaneous chiral symmetry breaking. 

Recent work on the description of bound states in these "mass-deformed" IR-conformal the- 
ories has clarified the situation, providing a framework for extrapolation and giving access to the 
mass anomalous dimension y* as a result [27, 59, 60, 61, 62, 63, 64]. The basic picture is as fol- 
lows: assume that we are simulating a theory with an infrared fixed point, and that the lattice cutoff 
A = has been tuned so that physics on scales below A is governed by the fixed point (in terms 
of the gauge coupling, o;(A) « a*). We now input a small fermion mass, m = m(A). The renormal- 
ized mass will then run, with its scale dependence determined by the mass anomalous dimension 
7*, 

m(jU) =m(A) . (4.1) 

As m increases, eventually it will reach a scale M such that m{M) = M. At energy scales below M, 
the now heavy fermions will decouple from the theory, and the remaining pure-gauge theory will 
cause the coupling to run and trigger confinement. Assuming that this happens over a relatively 
small range of energy scales, the "induced" confinement scale is therefore M, which is now a 
function of m: 

^ = „ji/(i+7*)^ (4.2) 

where I have now set a = 1 to focus on the mass dependence. A sketch of the overall situation is 
depicted in fig. 3. 

As in a QCD-like theory, the masses of bound states will be roughly set by the confinement 
scale; the main difference is that now the confinement scale itself vanishes as m — 0, with power- 
law dependence. All bound-state masses (including the pions, which are no longer Goldstone 
modes due to the lack of spontaneous symmetry breaking) will thus scale universally as 

Mx=cV/(i+^*). (4.3) 



Decay constants associated with these bound states are expected to scale in the same way [60, 61], 
although simulation results for have been observed to be qualitatively different. The chiral 
condensate is a special case, which is uniquely sensitive to ultraviolet scales and thus has several 
additional contributions; the leading quadratic divergence (^ff^f) ~ m/A^ is present in QCD calcu- 
lations as well. The next leading contribution arises from renormalization-group running between 
the UV and IR scales [63, 65, 62]: 

(VAVa) =Acm + Bcm^'^^-'^*'^l^^+'^*^^ + ... (4.4) 
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Figure 3: Sketch of the renormalization-group evolution in a mass-deformed theory with an IR fixed point. 
When a small fermion mass is introduced at UV cutoff A which is on the "plateau" a{A) « a*, the mass 
will run with the anomalous dimension at the fixed point: m{pL) ~ m{A){A/ jj.)"^ . At scale M such that 
m{M) = M, the fermions screen out of the theory and confinement is induced. The deformed theory has the 
usual spectrum of bound states, but the mass of any state scales as m'/('+''*' and vanishes in the limit m — > 0. 



Additional terms for this expansion are hypothesized in [63]. The precise behavior of the chiral 
condensate as a function of m remains a topic of some debate. The authors of [66] have studied 
the condensate in this context using the Schwinger-Dyson equation with the "ladder" truncation, 
giving some insight into the expected behavior and finite-size scaling. Recently, Patella has studied 
the condensate through the spectral density of the Dirac operator, finding that a relatively precise 
determination of y* can be made [67]. 

These "hyperscaling" relations have been applied to spectrum results for the SU(2) Nf = 2 
adjoint theory, with reasonable success [68, 69, 70, 71, 72, 73, 63]. Results for the mass anomalous 
dimension 7* vary depending on the method used, but are consistently much smaller than 1. With 
no spontaneous chiral symmetry breaking and a small 7*, this theory is of limited use for tech- 
nicolor model-building, but it remains a valuable proving ground for lattice methods designed to 
study infrared conformality. An interesting qualitative feature observed by [71] is that the glueball 
states in this theory are actually lighter than fermionic bound states, including the pseudoscalar me- 
son. Such an ordering was predicted as a possible signature of mass-deformed conformal theories 
simulated on the lattice [74]. 

The Nf = 12 theory with Nc = 3 and fundamental fermions has seen a great deal of interest, 
with a large number of independent groups carrying out simulations; a list of references is given 
in table 2. Although an eai^ly study of the SF running coupling indicated a relatively weak infrared 
fixed point [38, 39], there remains an ongoing debate in the community about whether Nf = 12 lies 
inside the conformal window. The set of groups studying the RG flow directly have generally found 
results consistent with an infrai^ed fixed point, while results for the spectrum of bound states are 
more contentious. As a particular example, simulation results from Fodor et al. [62, 75] have been 
fit by the authors using a simple polynomial ansatz and the mass-deformed functional forms shown 
here, with the conclusion that the IR-conformal fit is disfavored. However, the authors of [63, 64] 
argue for the opposite conclusion based on slightly different approaches to the mass-deformed fits. 
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Further and more rigorous development of the framework for dealing with mass-deformed CFTs 
will be required in order to resolve the disagreement. 

All of these discussions are further complicated by the observation of a nontrivial UV fixed 
point in the theory by [76, 77], which can obstruct the ability to extract continuum physics if 
simulations aie caiTied out too close to the alternate fixed point in bare parameter space. It is worth 
noting that all of the simulations carried out at Nf = 12 thus far have used staggered fermions. 
Certainly staggered fermions are an obvious choice for this theory, as no "rooting" is needed to 
thin out the number of continuum flavors as in application to QCD. However, some unusual results 
have been observed at strong coupling, in particular the breaking of single-site shift symmetry, 
manifesting as a difference between the plaquette as measured on adjacent lattice sites [78]. Many 
of the existing simulation results have been careful to avoid systematic effects induced by these 
strong-coupling artifacts, but the presence of such complications should be kept in mind and each 
study evaluated on a case-by-case basis. It may be worthwhile for lattice groups considering new 
studies of Nf = 12 to consider the use of Wilson or chiral fermions, to lend further clarity to the 
situation. 

Finally, I turn to the 3-color fundamental Nf = 10 theory. There has been much less work here 
than atNf = 12, but some preliminary results are beginning to appear I will discuss the results here 
in their own right, but note that these simulations cannot be completely separated from the debate 
at Nf =12. Any evidence that Nf = 10 is infrared conformal necessarily implies that Nf = 12 is 
within the window as well, and vice-versa if Nf = 12 is confining and chirally broken. 

The first results at Nf = 10 were obtained by [79, 80]. Their approach is to compute the 
Schrodinger Functional running coupling, using unimproved Wilson fermions and the plaquette 
gauge action. They follow the conventional approach of measuring the SF coupling g^{gQ,L/a) 
for many bai^e couplings and volumes, and then use step-scaling to reconstruct a discrete j3- 
function. This discrete jS -function is found to be consistent with zero for couplings between roughly 
10/3 < < 9.3. More recent results, which have been presented but not yet published, include 
a measurement of the mass anomalous dimension, pointing to a relatively large value between 
roughly 0.5<7*< 1 [81]. 

Very recently, the LSD collaboration has released spectrum results for the Nf = 10 theory 
[82]. The results are found to be well-described by a mass-deformed conformal fit, which yields 
a global best-fit value for the anomalous dimension consistent with 7* ~ 1, in agreement with 
the Schrodinger Functional result. However, only results for a single lattice volume 32^ x 64 are 
presently available, and finite-volume corrections (including topological-charge effects) remain an 
outstanding issue. 

If these preliminary hints of approximate scale invariance with 7* ~ 1 are found to be robust, 
then Nf = 10 will become a theory of great interest. In the context of technicolor models, it would 
(if outside the conformal window) be the first proven example of a "walking technicolor" theory, 
which could be used as the foundation for building a concrete model without issues with precision 
constraints. Even if Nf = 10 lies within the conformal window, it can be used to build a technicolor 
model by giving some of the fermions explicit masses, in order to induce confinement for the 
remaining fermions [27, 28]. In addition, if it is a theory with a strongly-coupled IR fixed point, 
the properties of A'y = 10 may be of interest from the perspective of AdS/CFT duality. In particular, 
a strongly-coupled CFT with a supergravity dual should exhibit a large splitting in the anomalous 
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dimensions of operators with spin greater than two, which correspond to 'stringy' excitations of the 
dual theory, from the lower-spin operators [83]. 

5. Outlook and conclusions 

The task of simulating many-fermion theories provides some novel challenges compared to 
dealing with QCD. Even far from the edge of the conformal window, the addition of fermions 
slows down the running of the gauge coupling, so that much larger changes in physical scale are 
needed to explore the same range of coupling strengths. Based on this we might expect that finite- 
volume corrections become much more problematic in many-fermion gauge theories. Indeed, the 
usual rule of thumb Mj^L > 4 for QCD simulations to have small finite-volume effects is found to 
be badly violated in an approximately scale-invariant system [84] . 

In addition to the requirement of large volumes to obtain accurate results, the presence of many 
fermion species further adds to the cost of HMC evolution. On top from a trivial scaling with Nf of 
the number of matrix inversions needed to compute the fermion force, the magnitude of the force 
also increases with Nf, so that the overall cost to maintaining fixed acceptance goes as Njr . All of 
these issues are compounded by the appearance of factors of Nf in the chiral perturbation theory 
expansion as discussed, which necessitates moving to lighter fermion masses (and thus even larger 
volumes.) 

The state of the field of many-fermion lattice simulations is similar in many ways to that of 
QCD simulations many years ago. Our understanding of these theories is growing rapidly thanks 
to the combined efforts of many groups, but precise study of near-conformal gauge theories will 
hkely require both increases in computing power and innovations in algorithms and techniques. In 
the meantime, answers from the LHC may help us to narrow down the parameter space to a few 
regions of interest. Even if the physics behind electroweak symmetry breaking is found to be a 
strongly-coupled gauge theory that has yet to be simulated, the work carried out thus far on models 
within and near the conformal window has been crucial in developing our ability to understand and 
work with such theories. 

Finally, although 1 have not discussed the topic here, it should be noted that there are several 
important ways that QCD simulations can contribute to searches for BSM physics. For decay pro- 
cesses which involve mesons and baryons, whether the decay is mediated by the standard model 
or by new physics, non-perturbative input is required in order to compute the expected rate. In 
flavor physics, where BSM theories generically predict a number of effects, uncertainty in QCD 
contributions remains a significant source of error; for a recent review, see [85], and [86] in these 
proceedings. In some cases, the interaction of BSM physics can require the computation of addi- 
tional operators that do not appear in standard model processes, including scalar and tensor matrix 
elements [87, 88] and right-handed currents [89, 90]. 

1 would like to thank the organizers for putting together a very successful conference, and for 
inviting me to review this topic. 1 would also like to thank Y. Aoki, C. Bouchard, M. Buchoff, 
T. DeGrand, G. Fleming, A. Hasenfratz, E. Itou, T. Karavirta, H.-W. Lin, M. Lin, H. Ohki, E. 
Pallante, A. Patella, D. Schaich, and G. Voronov for providing me with their results prior to my 
talk, and for many interesting discussions before and during the conference. Fermilab is operated 
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Table 2: Collection of references to simulation results aimed at studying the conformal transition in various 
classes of S\J{Nc) gauge theories. References are in roughly chronological order, except that subsequent 
results from a single lattice group are cited consecutively. The "other" column encompasses calculations 
of observables which are not part of the low-lying spectrum and are not used to determine RG flows, for 
example the S'-parameter. Results on Ns from these references are summarized in fig. 2. 
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